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Abstract 

Background: Transforming growth factor-beta (TGF-pi) gene has been regarded as an important mechanism in 
angiogenesis, endothelial cell proliferation, adhesion,and the deposition of extracellular matrix. The TGF-pi gene may be 
involved in the development of diabetic retinopathy (DR) through disrupting angiogenesis. However, studies investigating 
the relationship between — 509C/T and +869T/C(L10P) polymorphisms and DR yielded contradictory and inconclusive 
outcomes. In order to realize these ambiguous findings, a meta-analysis was performed to assess the association between 
the TGF-pi gene polymorphisms and susceptibility to DR. 

Methodology/Principal Findings: VJe conducted a search of all English reports on studies for the association between the 
TGF-pi gene polymorphisms and susceptibility to DR using IVledline, the Cochrane Library, EMbase, Web of Science, Google 
(scholar), and all Chinese reports were identified manually and on-line using CBMDisc, Chongqing VIP database, and CNKI 
database. The strict selection criteria and exclusion criteria were determined, and odds ratios (ORs) with 95% confidence 
intervals (CIs) were used to assess the strength of associations. The fixed or random effect model was selected based on the 
heterogeneity test among studies. Publication bias was estimated using Begg's funnel plots and Egger's regression test. 

Results: A total of three studies were included in the meta-analysis and all included studies analyzed patients with type 2 
diabetes. For -f869T/C(L1 OP) polymorphism, significant association was observed in an allele model (L versus P: OR =1.34, 
95%CI = 1.03-1.73) and the recessive model (LL versus LP-nPP: OR=1.70, 95%CI = 1.13-2.56). As regards -509C/T 
polymorphism, no obvious associations were found for all genetic models. 

Conclusions: This meta-analysis suggested that the H-869T/C{L1 OP) polymorphism in TGFpi gene would be a potential 
protect factor for DR. However, the -509CAr polymorphism is not associated with DR. 
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introduction 

It has been widely accepted that diabetic retinopathy (DR) is 
one of the foremost causes of bhndness in the working age 
population [1], characterized by angiogenesis in retina. However, 
the etiology of DR remains unknown and disease-modifying 
treatments are limited. In addition, since the involvement of 
cytokines in DR is hypothesized, there were many candidate genes 
approach in designing a case-control association study of single 
nucleotide polymorphisms (SNPs) including transforming growth 
factor-beta (TGF-pi) [2]. 

TGF-Pl has an important role in angiogenesis, endothelial cell 
proliferation, adhesion and the deposition of extracellular matrix 
[3,4]. The TGF-pi gene may be involved in the development of 



diabetic retinopathy (DR) through disrupting angiogenesis and 
blood retina barrier breakdown [5]. The highly polymorphic 
human TGF-pi gene is located on chromosome 19q 13.1-13.3 
[6]. There were some known TGF-Pl gene polymorphisms such 
as -988C/A, -800G/A, -509C/T and -h869T/C(L10P). 

Thus far, previous studies concerning association between 
— 509C/T polymorphism and risk of DR are limited and rather 
conflicting within type 2 diabetes [7,8,9]. The second TGF gene 
polymorphic exchange -l-869T/C(L10P) is not commonly analyzed 
as — 509C/T. Recendy, few findings were frequently performed 
on the effect of -H869T/C(L10P) polymorphism on DR, but the 
results were limited [7,8] . Considering that a single study may lack 
the power to provide a reliable conclusion, we performed a meta- 
analysis on these eligible studies, to investigate the precise 
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relationship between tiie TGF gene polymorphism and suscepti- 
bility to DR, which could have a much greater possibility of 
reaching reasonably strong conclusions. 

Methods 

Selection of eligible studies 

Two reviewers (Lei Liu and Jinghua Jiao) independendy 
scrutinized studies on the associations between TGFpi gene 
polymorphisms and DR. We searched Medline (Jan. 1st, 1946 to 
Oct. 31th, 2013), Embase Jan. 1st, 1950 to Oct. 31th, 2013), the 
Cochrane Library (up to 2013, issue 10), Chinese Biological 
Medicine Jan. 1st, 1978 to Oct. 31th, 2013), China National 
Knowledge Infrastructure [fan. 1st, 1979 to Oct. 31th, 2013), 
Wang Fang Data Jan. 1st, 1982 to Oct. 31th, 2013) and 
Chongqing VIP database Jan. 1st, 1982 to Oct. 31th, 2013) 
databases using the terms "transforming growth factor beta 1 or 
transforming growth factorpi or TGF beta 1 or TGFpi", 
"diabetes or diabetic complications" and "polymorphism, variant 
or mutation". We used the PubMed option "Related Articles" for 
each study to retrieve additional potentially relevant articles. 
Reference lists were checked and researchers contacted for 
additional literatures. Authors of publications ^\-ere contacted 
when results were unclear or when relevant data were not 
reported. The search was done without restriction on language, 
but we only included published articles written in English or 
Chinese. 

Selection criteria 

Studies were selectc'd if they met the following criteria: (1) 
association study in sporadic DR subjects; (2) there were available 
data for TGPpi gene mutations with risk of DR, using a case- 
control or cohort design; (3) the genotype distribution in the 
controls of all studies should be in agreement with Hardy- 
Weinberg equilibrium (HWE); and (4) in the case of multiple 
publications from the same study group, the most complete and 
recent results were used. 



the proportion of total variation contributed by between-study 
variation. If there was a statistical difference in terms of 
heterogeneity, the random effects model would be used to estimate 
the pooled ORs [11,12]. Otherwise, the pooled ORs were 
estimated using the fixed effects model [13]. Sensitivity analysis 
was carried out by deleting one single study each time to examine 
the influence of individual data set on the pooled ORs. The 
possible publication bias was estimated with funnel plots and 
Egger's test. An asymmetric plot suggests a possible publication 
bias and the P value of Egger's test less than 0.05 was considered 
representative of statistically significant pubhcation bias [14]. All 
statistical tests were performed with Comprehensive Meta-Analysis 
software version 2.0 (Biostat, Englewood Cliffs, I.N.J., USA) and 
RevMan version 5.0 (Review Manager, Copenhagen: The Nordic 
Cochrane Centre, The Cochrane Collaboration, 2010). All P 
values were two sided and a P value of smaller than 0.05 for any 
test was considered to be statistically significant. 

Results 

Study inclusion and characteristics 

Finally, a total of three studies fulfilling the inclusion criteria 
were identified [7,8,9]. The results of the search are shown in 
Figure 1. However, patients in these three included studies were 
type 2 diabetes. In one of these studies, cases were recruited from 3 
different countries, and thus each country in the literature was also 
considered separately for meta-analysis. Therefore, a total of five 
studies (three studies with 521 cases and 580 controls for — 509C/ 
T polymorphism and two studies with 268 cases and 340 controls 
for -l-869T/C(L10P) polymorphism) were included in the meta- 
analysis. The studies identified and their main characteristics are 
summarized in Table 1, Table 2 and Table 3. Genotype 
distribution of any polymorphism did not differ from Hardy- 
Weinberg equilibrium with in both groups (aU were greater than 
0.05). 



Exclusion criteria 

The exclusion criteria were defined as: 1) useless data reported, 
genotype number or frequency not included; 2) abstracts, reviews 
and animal studies; and 3) genotype distribution in the control 
population not consistent with HWE. 

Data extraction 

After excluding the overlap studies and including the additional 
ones, two investigators independently extracted data from each 
study with a standard fashion and entered into a common 
database. When discrepancies were appeared, all investigators 
were recruited to assess the data. The following information was 
collected: First author, year of publication, country, ethnicity, 
charactc-ristics, study design, sample sizes of patients and controls, 
genotype numbers, minor allele frequency (MAP) and numbers 
(MAN), P value for HWE. 

The review and analysis were guided to conduct by the 
PRISMA statement for preferred reporting of meta-analysis [10]. 

Statistical analysis 

Odds ratio (ORs) with 95% confidence intervals (CIs) for 
genotypes and alleles were used to assess the strength of association 
between TGPpi gene pol)'morphisms and DR. The ORs were 
performed for the allele contrasts, additive genetic model, as well 
as recessive genetic model and dominant genetic model, respec- 
tively. Heterogeneity was examined with f statistic interpreted as 



' 9'atides ideutified'though'database^ 



:::;;v:d^ 



3-aitil^ideotified'i;)ficated-i«movedt< 



4recordsexcludedbecauselackof- 
sufficient data, not English or 
Chinese language, not case-control- 
studyoroAortstudy.*' 



lTecoids-«3duded-because^lati(Hi- 
offl»HWE-4' 



Figure 1. Flow chart demonstrating those studies that were 
processed for inclusion in the meta-analysis. 

doi:1 0.1 371/journal.pone.00941 eO.gOOl 
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Quantitative data synthesis 

As shown in Table 4, the results showed the evidence of an 
association between the decreased risk of DR and +869T/C(L10P) 
polymorphism in an allele model (L versus P: OR =1.34, 
95%CI= 1.03-1.73) and the recessive model (LL versus LP+PP: 
OR= 1.70, 95%CI= 1.13-2.56). There was no significant diflFer- 
ence between-study heterogeneity (Figure 2). 

The association between — ,509C/T polymorphism and DR 
were also shown in Table 4. The results indicated no relationship 
of -509C/T polymorphism with DR risk (OR=1.18, 
95%CI = 0.96-1.45 for C versus T; OR= 1.23, 95%CI = 0.76- 
1.99 for CC versus TT; OR= 1.26, 95%CI = 0.97-1.63 for 
recessive model; OR=1.13; 95%CI = 0.72-1.79 for dominant 
model) (Figure 3). 

There was no significant publication bias according to Begg's 
and Egger's tests (Begg, p = 0.21; Egger, p = 0.71). 
% 

a. 

.i Sensitivity analysis 

^ In this study, we deleted every single study each time, in order to 

ffi examining the influence of the individual data set to the pooled 

~ ORs. According to sensitivity analytical results, we found that 

S there was no substantial modification of our estimates after 

>, exclusion of individual studies, indicating that the results were 

S stable (data not shown), 

o 

c 

^ Discussion 

■§ There are three TGF-(3 isoforms, TGFpi-3, which play 

^ important roles in regulating inflammation, cell growth and 

% differentiation [15]. Previous research has revealed that TGFpi 

^ ^ and TGFP2 were profibrotic, but TGFP3 might be aiitifibrotic 

[16]. Because there is little information on TGFP2 and TGFP3 
g ^ with DR, a pertinent analysis of isoforms other than TGppi in this 

meta-analysis has not been done. 

Previous research suggested that +869T/C(L10P) gene poly- 
morphism is associated with diabetic nephropathy in Chinese. The 
role of TGFpi in pathogenesis of diabetic nephropathy is 
pathologic tissue fibrosis leading to organ failure [1 7] . In addition, 
TGppi is believed to have an important role in the pathogenesis 
J £ of fibrotic diseases in the eye including proliferative vitreoretino- 

^ B. pathy [18]. The relationship between TGppi and the DR appears 

to be generation of exacerbation of angiogenesis and inhibiting the 
endothelial barrier function in the eye [19]. Therefore, the local 
activation of TGF-pi expression may play an important role in the 
Q "I development of the phase of DR. 

^ S Except for evaluating the expression of TGpp 1 in vivo or in 

=5 g vitro [20,21], there were few studies in DR patients that 

^ c susceptibility to angiogenesis may be correlated with the presence 

ffi ^ of particular alleles at the TGF-Pl locus. However, the results are 

~ -g- inconsistent and inconclusive due to limited sample size and 

c ^ different study populations. In previous study by Beranek M et al. 

c 3 [7], there is a marginally significant contrast by the — 509C/T 

g- 5 variant. But in the study by Paine SK et al. [9], — 509C/T was no 

i S g significandy associated with DR. 

I g ^ To the best of our knowledge, this is the first meta-analysis to 

i ^ o explore TGpp 1 polymorphism in development and progression of 

^ DR. To achieve a more reliable and comprehensive conclusion on 

Q ■§ 3- both variants, we used meta-analysis to assess the association 



o o 



u E 



^ -o a 

Q to _; 

2 2 E between the TGppi gene and DR risk on the basis of data from 

^ ^ three studies. We summarized a significant association for +869T/ 

S ^ C(LIOP) polymorphism, especially in an allele model (L versus P: 

|h5 or = 1.34, 95%CI= 1.03-1.73) and the recessive model (LL 

If '^ versus LP+PP: OR= 1.70, 95%CI= 1.13-2.56). Because the 

< L±J X! ' 

included studies were case-control studies, we did not perform a 
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Table 2. Distribution of — 509C/T genotypes and allele frequencies among DR of cases and controls, and P-values of HWE in cases 
and controls. 



First 

Auther Number Cases(N) Controls(N) HWE C Allele(N) 

(p*value) (p*value) 

case control for case for control control 

group group CC CT TT group CC CT TT group case group group 

Beranel< 73 172 44 24 5 0.49 79 80 13 0.23 112 238 

M [7] 

Buraczynsl<a 195 168 79 85 31 0.31 40 44 16 0.51 243 124 

M [8] 

Paine SK 253 240 172 70 11 0.27 153 73 14 0.19 414 379 

[9] 



Abbreviations: DR: diabetic retinopathy; HWE: Hardy-Weinberg equilibrium. 
doi:l 0.1 371 /journal.pone.00941 60.t002 

subgroup analysis by study design. However, the results indicated 
no relationship of — 509C/T polymorphism with DR risk 
(OR=1.18, 95%CI = 0.96-1.45 for C versus T; OR=1.23, 
95%CI = 0.76-1.99 for CC versus TT; OR =1.26, 
95%CI = 0.97-1.63 for recessive model; OR =1.13; 
95%CI = 0.72-1.79 for dominant model). The results from our 
study seem to indicate that the — 509C/T polymorphism may be 
not associated with DR risk. Because heterogeneity was found 
among the studies, we employed random-effect model. Then a 
sensitivity analysis was carried out by removing one study for each 
time and re-running the model to determine the effect on the 
overall estimate. The estimates changed quite little, indicating that 
the results were stable. Owing to the limited number of included 
studies, we did not perform a publication bias test. 

In addition, some other TGFpi genes including — 988C/A, 
-800G/A, and +915G/C(R25P) were studied in the relationship 
with DR. The — 988A allele was not represented in the study and 
no signicant difiFerences between groups were found for the 
-800G/A. However, +915G/C(R25P) polymorphisms in the 
TGF-Pl gene could be regarded as a strong genetic risk factor for 
DR [7]. As only one study was reported for in these candidate 
genes, we could not use meta-analysis to analyze the relationship 
between these genes and DR. 

In this meta-analysis, all included studies analyzed patients with 
type 2 diabetes. The etiology and mechanisms of type 1 and 2 
diabetes are different, so there may be diversity in the association 
between TGFpi gene polymorphism and DR with type 1 or 2 
diabetes. However, the relationship between TGFpi gene 
polymorphism and DR within type 1 diabetes is not demonstrated. 
More efforts are encouraged to explore this association. 



This meta-analysis has pooled all the available results from the 
case-control studies, which has significantly increased the statistical 
power. However, the results of the present meta-analysis should 
also be interpreted within the context of its limitations. First, DR is 
a multi-factorial disease from complex interactions between 
environmental exposures and genes factors. In this meta-analysis, 
we had insufficient data to perform an evaluation of such 
interactions for the independent role of TGFpi polymorphisms 
in DR development. Second, the major limitation of our study is 
the relatively limited available studies included in the meta- 
analysis. Thus, investigations involving large size of different races 
are necessary for a more reliable evaluation on their associations. 
Third, our meta-analysis is based on unadjusted estimates because 
of a lack of original data. Forth, although every effort was made to 
ascertain all appropriate publications, it is likely that some were 
missed or displayed erroneously. In addition, we did not consider 
studies published in languages other than English/Chinese or data 
presented in abstracted form. 

In spite of these limitations, for DR candidate genes, we believe 
that the positive locus identified in our systematic meta-analyses 
that warrant follow-up with high priority with a prospective trial. 

In conclusion, this meta-analysis suggested that the +869T/ 
C(LIOP) polymorphisms in the TGF-pi gene could be regarded as 
a strong genetic risk factor for DR. However, the — 509T/C 
polymorphism is not associated with DR risk. At the same time, 
this result should be interpreted cautiously. To better understand 
the potential mechanism for DR in humans, in the future, large 
well-designed epidemiological studies in the susceptibility of DR 
evidence are needed to confirm this association. It also wiU be 



Table 3. Distribution of +869T/C{L10P)(L10P) genotypes and allele frequencies among DR of cases and controls, and P-values of 
HWE in cases and controls. 



First Auther Number Cases(N) ^WE Controls{N) ^WE ^ Allele(N) 



case 


control 






(p*value) 
for case 








(p*value) 
for control 


control 


group 


group LL 


LP 


PP 


group 


LL 


LP 


PP 


group 


case group group 


Beranel< IVl [7] 73 


172 38 


27 


8 


0.35 


47 


94 


31 


0.18 


103 188 


Buraczynsl<a M 195 
[8] 


168 43 


86 


66 


0.14 


22 


44 


34 


0.28 


172 88 



Abbreviations: DR: diabetic retinopathy; HWE: Hardy-Weinberg equilibrium. 
doi:l 0.1 371 /journal.pone.00941 60.t003 
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Figure 2. Forest plot of the association between DR and the TGppi +869T/C(L10P) mutation (L vs P). Figure 26. Forest plot of the 
association between DR and the TGFpi +869T/C(L1 OP) mutation (LL vs PP). Figure 2C. Forest plot of the association between DR and the TGF|31 
+869T/C(L10P) mutation (LL vs LP+PP). Figure 2D. Forest plot of the association between DR and the TGFpi +869T/C(L1 OP) mutation (LL+LP vs PP). 
Figure 3A. Forest plot of the association between DR and the TGFpi — 509C/T mutation (C vs T). Figure SB. Forest plot of the association between DR 
and theTGFpi -509C/T mutation (CC vs TT). Figure 3C. Forest plot of the association between DR and the TGFpi -509C/T mutation (CC vs CT+TT). 
Figure 3D. Forest plot of the association between DR and the TGFpi — 509C/T mutation (CC+CT vs TT). 
doi:1 0.1 371 /journal.pone.00941 60.g002 



Table 4. Summary ORs and 95%CI of the association between polymorpliisms in the TGF (51 genes (-509CAr and +869T/ 
C(L10P)(L10P)} and DR risl<. 





Gene 


Allele contras 




Additive model 




Recessive model 


Dominant model 




OR(95%Cl) 


P 


OR(95%CI) 


P 


OR(95%CI) 


P 


OR(95%Cl) P 


+869T/C(L10P) 


1.34(1.03-1/73) 


0.03 


1.55(0.93-2.60) 


0.1 


1.70(1.13-2.56) 


0.01 


1.19(0.78-1.83) 0.42 


-509C/T 


1.18(0.96-1.45) 


0.11 


1.23(0.76-1.99) 


0.4 


1.26(0.97-1.63) 


0.09 


1.13(0.72-1.79) 0.6 



Abbreviations: DR: diabetic retinopathy; ORs: odds ratios; CI: confidence intervals. 
doi:10.1371/journal.pone.0094160.t004 
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Figure 3. A. Forest plot of tiie association between DR and tiie TGppi -509C/T mutation (C vs T). Figure 3B. Forest plot of the association 
between DR and the TGFpl -509C/T mutation (CC vs TT). Figure 3C. Forest plot of the association between DR and the TGF|31 -509CAr mutation (CC 
vs CT+TT). Figure 3D. Forest plot of the association between DR and the TGFpi -509C/T mutation (CC+CT vs TT). 
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necessary to combine genetic factors and other environmental risk 

factors. 
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